Nucleation of solitary wave complexes in two-component mixture Bose-Einstein 

condensates 
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Two-component mixtures of Bose-Einstein condensates are shown to support solitary wave com- 
plexes that move with a constant velocity preserving their form. I obtain the families of such solitary 
wave complexes in two-dimensional two-component mixture condensates. These solutions are classi- 
fied according to the structure of the wavefunction in each component. I show that these complexes 
nucleate from the surface of the disk when it moves supercritically, therefore, suggesting a mecha- 
nism by which these waves can be obtained in condensates by a moving laser beam. The condition 
for such a nucleation is derived analytically. The fiow for supercritical disk velocities is computed 
numerically. The process of a boundary layer separation with emission of either vortex pairs in 
each component or a vortex pair in one component and a "slaved wave" in the other component is 
elucidated. 

PACS numbers: 03.75.Lm, 05.45.-a, 67.40.Vs, 67.57.De 
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The systems of the coupled nonlinear Schrodinger 
(NLS) equations are fundamental and universal systems 
that have been used to describe motions in conserva- 
tive systems of vifcakly nonlinear dispersive waves in 
continuum mechanics, plasma physics, nonlinear optics 
and condensed matter. Recent experimental advances 
in multi-component Bose-Einstein condensates (BECs) 
stimulated the interest in solitary wave solutions of these 
equations and their production as they define the possible 
excitation states that multi-component BECs can sup- 
port. Multi-component condensates have been formed 
by simultaneous trappiiig and cooling of atoms in distinct 
spin or hyperfine levels . There is also a hope to obtain 
condensed mixtures of different atomic species ^] . In this 
respect a lot of attention has been paid to discovering the 
topological solitons and defects such as domain walls 
and skyrmions (vortons) 0] in the regime of phase sep- 
aration. In the regime where phases do not separate, 
the complete families of three-dimensional solitary wave 
solutions were obtained in . There solutions were clas- 
sified according to the structure of the wavefunction of 
each component. Four basic types were noted: (1) vor- 
tex rings of various radii in each component, (2) a vortex 
ring in one component coupled to a rarefaction solitary 
wave of the other component, (3) two coupled rarefaction 
waves, (4) either a vortex ring or a rarefaction pulse cou- 
pled to a localised disturbance of a very low momentum. 
The continuous families of such waves were shown in the 
momentum-energy plane for various values of the inter- 
action strengths and the relative differences between the 
chemical potentials. 

Two-dimensional solutions of the coupled NLS equa- 
tions are relevant in the view of the experimental setting 
for studies of multi-component BECs, for instance, two- 
dimensional vortex complexes were shown to nucleate in 
rotating two-component condensates giving rich equi- 
librium structures such as triangular, square and double- 



core lattices and vortex sheets. The goal of this letter 
is to obtain complete families of solitary wave solutions 
in two-dimensional two-component mixture BECs and to 
find the conditions for their creation in bulk condensates. 

For two components, described by the wave functions 
ipi and ip2 , with A^i and N2 particles respectively, the sys- 
tem of the Gross-Pitaevskii (CP) equations (the coupled 
NLS equations) on the wave functions of the condensate 
components become 
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where 



is the mass of the atom of the ith conden- 



sate, and the coupling constants Vij are proportional 
to scattering lengths a.y via Vij = ^ttTi^ a.ij / rriij , where 
uriij — rriimj / {rrii + ruj) is the reduced mass. The inter- 
action conserves the number of atoms of the two species, 
so / iV'iP dxdy — Ni. To study the equilibrium properties 
the energy functional has to be minimized subject to con- 
straint on conservation of particles leading to introduc- 
tion of two chemical potentials ni = VnUi -f Vi2n2, /i2 = 
Vi2ni + V22«2, where rii = \^pi\'^ is the number density 
in equilibrium. The dispersion relation between the fre- 
quency uj and the wave number k of the linear perturba- 
tions ((X exp[zk • X — iut]) around homogeneous states is 
obtained as 
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where uJ~{k) = cfk^ + fi k /Amj coincides with a one- 
component Bogoliubov spectrum with the customary de- 



fined sound velocity 



niVii/mi and — 2^2^ 



where — nin2Vi2/mim2- The system is dynami- 
cally stable if the spectrum J^l is real and positive which 
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implies that V11V22 > for stability. The acoustic 
branches of Eq. (jJl are uj± ^ c±k with the correspond- 
ing sound velocities 2cj_ = + c| ± {c\ — c^Y + 4c|2- 
The solitary waves I seek below are all subsonic, so their 
velocity U is less than c_. 

Solitary waves. I shall restrict the parameter space by 
letting nil = 7712 = m, Vn = V22, a = Vi2/V"ii. To 
find axisymmetric solitary wave solutions moving with 
the velocity U in the positive a;— direction, I solve 

2i[/^ = V^Vi + (1 - - a\4,2?)^i (4) 
2iJ7^ = V2V'2 + (1 - al^iP - |^2p - A2)^2, (5) 

V'l^V'lcx), -02 ^"0200, as |x| ^ 00, 



FIG. 1: (colour online) The dispersion curves of several fami- 
lies of the solitary wave solutions of Eqs. (0]]SJ . The solid lines 
show three families of solutions with a = 0.5 and = 0.1. 
The numbers next to the dots give the velocity of the soli- 
tary wave. The top (black) branch corresponds to VP- VP 
complexes. The middle (green) branch shows p ys £ for VP- 
SW complexes and the bottom (red) branch is the dispersion 
curve of SW-VP complexes. The dashed (blue) line across the 
solid black branch shows the VP- VP complexes for U = 0.2, 
= 0.1 as the intercoupling parameter a increases in in- 
crements of 0.1 from 0.1 (top point) to 0.7 (bottom point). 
The light grey (magenta) dashed line shows the VP- VP com- 
plexes for U = 0.2 and a = 0.5 with the asymmetry param- 
eter A^ taking values 0.05 (top point), 0.1,0.2,0.3 (bottom 
point). The insets show the plots of 2; = |7/;i(a;, j/)p (top) 
and z — \tlj2{x,y)\'^ (bottom) for the VP-SW complex with 
U = 0.2, a = 0.5, A^ = 0.1. 



where a dimensionless form of Eqs. is used, such 

that the distances are measured in units of the corre- 
lation (healing) length ^ = Ti/^/Sm/Tj", the frequencies 
are measured in units 2fii/h and the absolute values 
of the fields and |'02p are measured in units of 

particle density fJ-i/Vn. Also present in Eqs. (^HS is 
the measure of asymmetry between chemical potentials 
= {fii — /i2)/Mi (where we assume that fii > /i2). 
The condition of dynamical stability under which two 
components do not separate is < 1. The values of 
the wave-functions of the solitary waves at infinity in 
Eqs. (gHHIl are given by ip^oo = (1 - a - A2)/(l - a^) 
and i/'ioo = 1 ~ Q^'^/'ioo and the critical speed of sound is 

that c_ > implies A^ < 1 — a. 

To find each solitary wave complex, I will fix a,A2, 
and U, solve Eqs. (@HH1) by Newton-Raphson iteration 
procedure described in 5] and calculate the momenta 
Pi = 11 /[(0i" ~ 0ioo)V0i - (0., - 0ioo)V0*] dxdy and 
energy 

2 

^ - I /E{iv0,p + ^(0L-i0.pnrf^rfy 

i—l 

+ f / t[{^L-\i^^\')dxdy (6) 
•' i=i 

per unit length in z-direction. Similarly to the three- 
dimensional case we can show that U = dS/J^Pi 
where the derivative is taken along the solitary wave se- 
quence. Also, we multiply Eqs. by xdij^l/dx and 
their complex conjugates by xdipi/dx, integrate by parts 
over all space, and compare the result with © which 
gives £ — /^|90i/9x| dxdy. Similarly, firstly we re- 
place 0i by 0i — 0'ioo in the first term of the first integral 
of Eq. (jnj and integrate by parts, secondly we multi- 
ply Eqs. by ydip* /dy and integrate by parts which 
gives us two more integral properties that will be used as 
checks of the numerical work 
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+«y'n(v^?oo-i0.i')'^^dy, (7) 

f = ^/(i-i0ip-"i02n 

x(30i2^ - 0ioo(0i + 01*) - |0in dxdy 

+ iy'(l-a|0i|2-|02|^) 

x(30|oo - 02oo(02 + 02*) - \i>^?)dxdy. (8) 

In the limit a — > 0, two components become uncoupled, 
in which case the solitary wave sequence for each compo- 
nent follows the dispersion curve of the one-component 
GP equation calculated in 0. The family of sohtary 
waves in 2D is represented by a pair of point vortices 
(VP) of opposite circulation if Uc < O.SGq. These vor- 
tices are separated by distance 2bi ^ U^^ for small U. 
As the velocity increases, the wave loses its vorticity and 
becomes a rarefaction pulse (RP). As U ^ Ci both en- 
ergy and momentum per unit length approach zero and 
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the solutions asymptotically approach the 2D rational so- 
lution of Kadomtsev-Petviashvili Type I equation. The 
sequence merges tangentially with the phonon branch of 
the dispersion curve in each of the uncoupled compo- 
nents. For a ^ ci ^ C2, so different components be- 
come RP at different critical values of U and a variety 
of complexes becomes possible. Table 1 gives an exam- 
ple of various transitions from one complex to another as 
the velocity U increases in the system with a = 0.05 and 
= 0.1. 

Table 1. The velocity, U, energy, momenta, pi, and half- 
separations between centres of the point vortices, bi, of the solitary 
wave solutions of Eqs. |4][^ with a = 0.05 and = 0.1. The 
sequence terminates at (7 = c_ fs; 0.646. 

U £ pi p2 bi 62 complex 
0.40 14.7 13.8 12.1 0.915 0.498 VP- VP 
0.43 13.7 12.5 10.9 0.184 - VP-RP 
0.45 13.0 11.7 10.2 - - RP-RP 
0.5 11.4 9.90 8.46 - - RP-RP 
0.6 7.68 6.29 5.36 - - RP-RP 

Next I will consider the cases of more intermediate val- 
ues of intercoupling interaction strength. In addition to 
the VP- VP, VP-RP, and RP-RP complexes, there are 
new classes of solitary waves that have no analog in one- 
component condensates. In these complexes the distur- 
bance of a very low impulse in one condensate (called 
"slaved wave" (SW)) is dragged by either VP or RP 
structure of the other component. The density of SW is 
maximal where the density of either VP or RP is minimal 
and vise versa. For fixed values of a and the system 
has three families of solitary wave complexes: VP(RP)- 
VP(RP), SW-VP(RP) and VP(RP)-SW as Fig.l illus- 
trates. Also, Fig^ shows the dispersion curves of several 
other families of the solitary wave solutions in the system 
when two out of three parameters (a, A^ and U) are kept 
fixed. 

Nucleation. In a pioneering paper Frisch et al. Q used 
a direct numerical simulation of the one-component GP 
equation to show that the superflow around a disk re- 
leases vortices from the perimeter of the disk creating a 
net drag force beyond a critical velocity. The criteria for 
vortex nucleation was related to the transonic transition, 
namely, the vortices are created when the local speed of 
sound is reached somewhere in the mainstream. The ar- 
gument supporting this conclusion is based on the obser- 
vation that the hydrodynamical form of the steady state 
one component GP equation away from the disk bound- 
ary changes its type from being elliptic to hyperbolic if 
the local speed of sound is exceeded. One would expect 
that somewhat similar scenario should exist for vortex 
nucleation in multi-component condensates. But in two- 
component condensates there are several sound speeds: 
ci,C2,ci2,c_ and 0+. So which of these values leads to 



vortex nucleation? The argument based on the criterion 
of the dynamical stability that the vortices should nucle- 
ate as soon as the flow reaches the local c_ value can be 
refuted by the reference to a one-component GP equa- 
tion with a nonlocal potential that allows for the roton 
minimum in the dispersion curve 0] and therefore has 
two critical velocities: speed of sound c and the Lan- 
dau critical velocity vl, vl < c. We have shown in 
that the nucleation of vortices is related to c and not to 
Vl- Another difficulty is that in multi-component con- 
densates each component has its own velocity u^. Will 
vortices nucleate when the velocity of just one compo- 
nent reaches the criticality or is there a more intrinsic re- 
lationship between the velocities of the components and 
critical velocity? Finally, as we discovered above there 
exist more than one vortex complex, so which complex 
nucleates at the criticality? Next I answer these questions 
using a simple analytical argument and direct numerical 
simulations. 

By using the Madelung transformations V'i — 
Ri exp[iS'i] in Eqs. and separating the real and 

imaginary parts, one gets the following hydrodynami- 
cal equations for the number density Ui — i?f and the 
phase 4>i = fiSi/m for the superflow with 0,^ — UqoX as 

+ — > 00 

^ + V • (n,V0,) = (9) 
d(j), 1 2 1 2 , / N 

+ A^Z^, (10) 

where j = 2 if i = 1 and j — 1 if i — 2. We con- 
sider a stationary flow and neglect the quantum pressure 
terms on the right hand side of Eq. H10|l due to our 
interest in mainstream flow. Since two components are 
coupled through their amplitudes, the velocity vectors of 
two components are parallel. We flx a point outside of 
the disk at which the components move with velocities 
Ui, introduce the local orthogonal coordinates such that 
the a;— axis is tangent to the flow and expand in the 
neighbourhood of this point as « UiX + e(f>i, where e 
is a small parameter. To the leading order Eqs. H9I10() 
become 

-^4^ixxxx ~t~ ^4^ixxyy ~t~ C(l)iyyyy — 0, (f f ) 

where A — d{niui) / dui x d{n2U2) / du2 ~ uiU2dni / du2 x 
dn2/dui, B = n2d{niui) / dui + nid{n2U2) / du2 and C = 
nin2 where now rii = ni{ui,U2). Note, that > AAC 
for all Ui, therefore, Eq. Hll|l is elliptic (in the sense that 
it has no real characteristics d(f>i/dy + Xd(f>i/dx = 0) 
if and only if A > 0. This gives the criterion for the 
vortex nucleation: the boundary layer separation with 
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nucleation of vortices takes place when 

d{niui) d{n2U2) dni dn2 

—5 H =WlW2T^T^, (12) 

oui au2 0U2 oui 

somewhere in the mainstream. In our special case of 
equal masses and intracomponent coupling parameters, 
condition (|12(l in dimensionless units becomes 

{ni-2u\){n2~2ul) ^ Aa^n\nl. (13) 

Note that if a = 0, then the criterion (|13|l says that the 
criticality occurs when the mainstream velocity reaches 
the local speed of sound 2uc = ric- If a 7^ and if 
at criticality ui ~ U2 = Uc, then Uc — ci-., where q_ 
is a local speed of sound defined by 4cf_ = rii + ^2 — 

(ni — n2)^ + 4Q!^nin2, but in general Eq. 1)13(1 can be 
satisfied when ui and U2 bracket q_. 

Direct numerical simulations of Eqs. (^1211, with —2\il)it 
added to the left hand sides and in the frame of reference 
in which the disk is stationary so U = u^o, show vortex 
complex nucleation in supercritical flow around the disk. 
This suggests that these complexes could be generated 
by a laser beam which moves supercritically in trapped 
condensates. At subcritical velocity {U < 0.225), the 
flows of the condensates are symmetric fore and aft of 
the direction of motion, and the disk experiences no drag. 
When the condition p3|l is satisfied, which happens first 
on the disk equator where the velocities are maximal, the 
condensates evade shocks through a boundary layer sep- 
aration. Fig. |2] shows the emission of various complexes 
for the disk of the radius 10 healing lengths that moves 
with supercritical velocity U = 0.28. The disk sheds 
SW-VP, VP-SW and VP- VP complexes in the order and 
frequency that depends on the value of the disk's ve- 
locity. These complexes move more slowly than the disk 
and form a vortex wave street that trails behind it, main- 
tained by other complexes that the disk sheds. As the 
velocity of the disk increases such a shedding becomes 
more and more irregular. Each complex is born at one 
particular latitude within the healing layer on the disk. 
As it breaks away into the mainstream, it at first con- 
tributes a flow that depresses the mainstream velocities 
on the disk below critical. For larger values of the disk 
velocity {U > 0.265), more energy is required for this 
depression and VP- VP complex that has larger energy 
than SW-VP complexes is born first. At low supercriti- 
cal velocities SW-VP complex is born first. As it moves 
further downstream however, its influence on the surface 
flow diminishes. The surface flow increases until it again 
reaches criticality, when a new complex is nucleated and 
the whole sequence is repeated. The vortex and slaved 
wave street trailing behind the disk creates a drag on the 
disk that decreases as the nearest complex moves down- 
stream, but which is refreshed when a new complex is 
born. The complexes downstream of the disk move with 
different velocities and interact among themselves. These 



FIG. 2: (colour online) The time snapshots of lipi]^ (top) and 
|i/'2p (bottom) of the solution of Eqs. with —2iipit added 

to the left hand sides with a = 0.5, =0.1 for the flow 
around a disk of radius 10 moving to the right with velocity 
U = 0.28. The solitary wave street is seen in the wake of 
the disk. The complexes were emitted in the order VP- VP, 
SW-VP, VP-SW, and SW-VP complex has just got emitted 
from the disk boundary. For this large supercritical velocity 
the VP- VP complex is first nucleated from the surface of the 
disk and the insets show this moment at an earlier time. On 
the main panels this complex is in the process of splitting 
into the VP-SW and SW-VP complexes. Only parts of the 
computational box are shown. 




interactions may lead to a transformation from one type 
of the complex to another; Fig. |2 shows the splitting 
of the VP- VP complex into SW-VP and VP-SW com- 
plexes. The mechanism in which solitary waves transfer 
energy from one to another was elucidated in ^3 for a 
one-component condensate. 
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